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ABSTRACT 

Cosmological parameters from WMAP 7 year data are re-analyzed by substituting a 
pixel-based likelihood estimator to the one delivered publicly by the WMAP team. 
Our pixel based estimator handles exactly intensity and polarization in a joint man- 
ner, allowing to use low-resolution maps and noise covariance matrices in T,Q,U at 
the same resolution, which in this work is A s ;d c = 16. We describe the features and 
the performances of the code implementing our pixel-based likelihood estimator. We 
perform a battery of tests on the application of our pixel based likelihood routine to 
WMAP publicly available low resolution foreground cleaned products, in combination 
with the WMAP high-£ likelihood, reporting the differences on cosmological parame- 
ters evaluated by the full WMAP likelihood public package. The credible central value 
for the cosmological parameters change below the 1 a level with respect to the evalu- 
ation by the full WMAP 7 year likelihood code, with the largest difference in a shift 
to smaller values of the scalar spectral index ris- 

Key words: Cosmology: cosmic microwave background - Physical data and processes: 
cosmological parameters. 



1 INTRODUCTION 

The anisotropy pattern of the cosmic microwave background 
(CMB) is a treasure for understanding the costituents of 
our Universe and how it evolved from the Big Bang. Un- 
der the assumption of isotropy and Gaussianity of CMB 
fluctuations, the power spectra of intensity and polariza- 
tion anisotropies include all the compressed information on 
our Universe through the determination of the cosmologi- 
cal parameters. There has been a tremendous improvement 
in the estimate of cosmological parameters driven by the 
increasingly better quality of CMB data, mainly due to 
the full sky observations in temperature and polarization 
by t he Wilkinson Microwave Anisotropy P robe (WMAP), 
(see lLarson et all i|2010h ; iKomatsu et all (|201Ch and ref- 
erences therein) and to the small angular scales mea sure- 
ment s by QUAD in polarization dBrown et al.l 2009|). by 



SPT jLueker et alJ[2ulp|:lKeisler et al.l201ll;lReichardt et all 



l2012r i and ACT l|Das et al.ll201ll : iDunklev et al.ll201ll ) in 

temperature. Planck will lead to a drastic improvement 
of CMB full sky maps in temperature and polarization, 
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leading to an eagerly expected improvement in cosmo- 
logical parameters with un certainties at the percent level 
i|Planck Collaboration! 120051 ). 

A joint likelihood analysis in temperature and polar- 
ization is one of the accepted methods in securing the sci- 
entific expectations of observational achievements in terms 
of cosmological parameters. Although the likelihood could 
be written exactly in the map domain under the Gaus- 
sian hypothesis, its computation is almost prohibitive al- 
ready at the resolution of 2 degrees, whereas cosmologi- 
cal information is encoded in the temperature and polar- 
ization power spectra up to the angular scales of the or- 
der of few arcminutes, where the Silk damping suppress the 
CMB primary anisotropy spectrum. It is now commonly ac- 
cepted to use an hybrid approach which combines a pixel 
approach at low resolution with an approximated likeli- 
hood bas ed on power spectrum estimates at high mult i- 
poles (seelBond. Jaffe and Knoxl <|2000l h IVerde et al l (|2003T ): 
lHamimeche and Lewis! |2008) for some of these approxima- 
tions). 

Since the three year release of the full polarization 
information, the WMAP team adopted such a hybrid 
sch eme approa c h, wh ic h has been suggested independentl y 
by lEfstathioul (l2004h : ISlosar. Seliak and Makarovl (|2004l h 
lO'Dwver et alj l|2004l ). At a first appearance of the three 
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year data, the WMAP team a dopted a pixel approach on 
HEALPIX jGorski et alj|2005h resolution N sidc = 8 tem- 
perature and polarization maps, and considered the high-f 
approximated likelihood to start at I = 13 in temperature 
and £ = 24 in polarization and temperature-polarization 
cross-correl ation for the deter mination of cosmological pa- 
rameters in lSpergel et al. | \200H ). The WMAP team treates 
separately tempe ratur e and polariza t ion a s explained in 
IPaee et al l l|2007l ) and iHinshaw et all <|2007h . by using the 
approximation that the noise in temperature is negligible. As 
a consequence, the WMAP likelihood code includes either 
(Q, U) and the temperature-polarization cross -correlation in 
the sa me sub-matrix. It was then shown by lEriksen et al.l 
< |2007t ) that by increasing the resolution of the temperature 
map to HEALPIX N s id c = 16 and therefore the multipole 
of transition to high-l 1 approximated likelihood in tempera- 
ture from £ = 12 to I = 30, the mean value for the scalar 
spectral index n a shifted to higher values by a 0.4 a. The 
asymmetric handling of the low-resolution temperature map 
at A^idc = 16 and polarization at A^idc = 8, became the final 
treatment of the three year data release. This low-^ likeli- 
hood aspect in the WMAP hybrid approach has not changed 
since the final release of the WMAP 3 year data to the cur- 
rent WMAP 7 year one. 

In this paper we wish to perform an alternative determi- 
nation of the cosmological parameters from WMAP 7 pub- 
lic data, substituting the WMAP low-<? likelihood approach 
with a pixel based likelihood code which treats T, Q, U at the 
same HEALPIX resolution AT s ide = 16 connected to the stan- 
dard WMAP high-^ package. In this analysis we therefore in- 
crease the resolution of polarization products digested by the 
pixel base likelihood f rom A^ide = 8 to A/si de = 16, in anal- 
ogy with what done bv lEriksen et alJ l|2007l ) for temperature 
only. The WMAP 7 year foreground cleaned (Q, U) maps, 
covariance matrices and masks at the resolution A s id e = 16 
are also publicly available at http://lambda.gsfc.nasa.gov: 
therefore, all data used in this paper are made available by 
the WMAP team. 

The paper is organized as follows. In Section II we 
briefly describe the WMAP hybrid approach to the likeli- 
hood, with particular care to the low multipole part. In Sec- 
tion III we describe our pixel approach, implemented in the 
BoPix code. We then present in Section IV the cosmological 
parameters obtained by using our alternative pixel approach 
in place of the WMAP one for a ACDM scenario. In Sec- 
tion V we extend our investigations to other cosmological 
models. In Section VI we draw our conclusions. 



2 A BRIEF DESCRIPTION OF THE WMAP 
HYBRID LIKELIHOOD ANALYSIS 

In the map domain, the likelihood as function of the cosmo- 
logical parameters {8} 



£(d|0) = 



2ttC 



1/2 



exp 



-id'c-d 



(1) 



where the data, d = s + n, is a CMB fully polarized map, 
considered as a vector combining T, Q and U maps, the sum 

1 The number of pixels in a map is given by Af p ; x = 12N? dc , i.e. 
768 for 7V sidc = 8 and 3072 for 7V sidc = 16. 



of signal s and noise n; the quantity C = S + N is the total 
covariance matrix, the sum of the CMB signal covariance 
matrix 8(6*), and the noise matrix N. The signal covariance 
is determined by the powe r spectra Cf Y , where X,Y are 
any o f T,E,B as given in iTegmark and de Oliveira-Costal 
(|200lf ). 

The WMAP low - l like lihood is described in the Ap- 
pendix of iPage et all (|2007f ) and we report here the essen- 
tials. The WMAP approach is based on the assumption to 
ignore the noise in temperature, which leads to a simplifica- 
tion of the likelihood, useful from the numerical computation 
perspective. By assuming that the noise in temperature is 
negligible at low multipoles, the WMAP approach consists 
in rewriting Eq. fTJ as: 



C(d\i 



exp ( 



exp[-|d t P (^p + 7vp)- 1 dp] 
\S P + Np\^ 2 



(2) 



where St is the temperature signal sub-matrix, the new po- 
larization data vector is dp — sp + rip, with sp = (Q, U) 
given by 



^ tp qTE 1 
Q ~ j y y (jTT ^ ] a l m i + ^Ylrn + -2Y lm ), 



u 



■ l P (jTE f 



(3) 
(4) 



with Sp (Np) is the signa l (noise) covarian ce matrix for the 
new polarization vector (|Page et al.l |2007| ). The noise co- 
variance matrix for (Q, U) equals the original one for (Q , U) 
when the noise in temperature is zero |Page et al.ll2007M . As 
temperature aj^, the full-sky internal linear combination 
(ILC) map is used (Hinshaw et al. 2007). 

According to IPage et all (|2007i ). Eq. Q and Eq. (0 
are mathematically equivalent when the temperature noise 
is ignored. With this assumption, the new form, Eq. p]l. 
allows the WMAP approach to factorize the likelihood 
of temperature and polarization, with the information in 
their cross-correlation, Cf E , retained in the polarization 
sub-matrix. As already mentioned in the introduction, 
temperature is considered at the HEALPIX resolution 
A^ide = 16, whereas polarization at N B id e = 8. The range 
of multipoles considered in the polarization sub-matrix is 
up to the Nyquist limit at N B id a ~ 8, i.e. ip = 23. Two 
computation options are ava ilable for the temperature 
likelihood, Gibbs sampling Jjewell. Levin and Andersonl 



2004 Wandclt, Lars on and La kshmhiaravaiianl 120041 : 
Eriksen et al] l2004h with a range of multipole considered 

up to It = 32 and pixel base evaluatio n, with It = 3 0. 

T he high-l likelihood , described in lLarson et al.l (|20ld ) 
and in I Verde et al.l l|2003l ). has been updated to beam/point 
sources uncertai nties through the various subsequen t 
WMAP releases ijHinshaw et al] 120071 : iNolta et all 120091 ). 
The high-t? TT likelihood takes into account multipoles from 
I = 31 (£ — 33) when connected with the pixel (Gibbs) like- 
lihood evaluation of the low resolution temperature data up 
to £ = 1200; the high-*! TE (and TB whe n used) likelihood 
takes into account multipoles from £ = 24 l|Page et al]|2007l ) 
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to I — 800. High-£ EE and BB data have not used so far in 
the various relases of the WMAP likelihood code. 



3 BOPIX 

BoPix computes the likelihood function in Eq. (fTJ for the pa- 
rameter space {8} which the Cf Y ({6}) depend on, without 
any approximation and with the same resolution in temper- 
ature and polarization. BoPix is a multithreaded OpenMP 
Fortran90 library which can be con nected to a sampler - to 
CosmoMC iLewis and Bridlel (|2002h in this work. 

The computation of the likelihood given in Eq. (1) re- 
quires an environment initialization, in which BoPix calcu- 
lates the geometrical functions dependent on the cosine of 
the angle between two pixels and reads the noise covariance 
matrix (C-binary format). 

BoPix then starts to compute the signal covariance ma- 
trix S for a given Cf Y {{9}) with a OpenMP routine with a 
high intrinsic level of parallel architecture, to which the noise 
covariance matrix N is summed. The full covariance matrix 
is then Cholesky decomposed. The computation of the de- 
terminant is obtained from the properties of the Cholesky 
decomposed matrix L: detC = (detL) 2 . The term C _1 d 
is computed as the solution for the variable x (vector with 
dimension 3A p i x ) of the equation Cx = d. 

The matrix manipulations are implemented on LA- 
PACK and BLAS mathematical libraries (as nag, essl, acml 
and mkl). There is an effort to improve the BoPix capabil- 
ities and performances (in terms of run time and memory) 
to make the direct likelihood evaluation at low resolution 
for cosmological parameters extraction as fast as possible, 
in particular by reducing the time spent for the Cholesky 
decomposition, and optimizing the combined scalability in 
memory and CPU time of this code; indeed, the resources 
required by BoPix are larger than those for the WMAP low- 
£ likelihood code since the polarization sector is treated at 
higher resolution. At present, BoPix can handle maps and 
full noise covariances up to HEALPIX Aside = 32 resolu- 
tion. On IBM Power6 (4.2GHz) architecture, available at 
CINECA (http://www.cineca.it ), with 64 threads on 64 log- 
ical CPUs (32 cores) BoPix can calculate the likelihood in 
about 0.3 seconds at A s id c = 16, and in about 15 seconds at 
Aside = 32. At A s idc = 16 on the same IBM Power6, a good 
trade off between computation time and memory required is 
obtained for 2 sec with 8 cores. More details about perfor- 
mances an d compari s on am ong different platforms will be 
provided in lDe Rosal (|2012l V 



4 DATA SET FOR BOPIX 

We use the temperature ILC map smoothed at 9.812 5 de- 
grees and reconstructed at HealPix (|Gorski et al.| [20051 res- 
olution A s idc = 16, the foreground cleaned (unsmoothed) 
low resolution maps and the noise covariance matrix 
in (Q, U) publicly available at the LAMBDA website 



weigthing accordingly to the WMAP team |jarosik et al.l 
12007ft 



http:/ /lambda.gsfc. nasa.gov/ for the frequency channels Ka 
(23 GHz), Q (41G Hz) andV (61 GHz) as considered by 



dpoi = c po i(c K ^dA- a + Cq dQ + c v d 



(•>) 



where di, Ci are the polarization maps and covariances, re- 
spectively (for i=Ka, Q and V). The total noise covariance 
matrix is therefore: 



~pol 



-i . — i , -l 

C Ka T C Q + C V 



(6) 



This polarization data set has been extended to temperature 
consideri ng the ILC map with an extra noise term, as sug- 
gested in lDunklev et all |200|). We have therefore added to 
the temperature map a random noise realization with vari- 
ance of <Ty T = IfiK 2 and consistently, the noise covariance 
matrix for TT is taken to be diagonal with variance equal 
to l^iA" 2 . The total noise covariance N for WMAP 7 yr data 
is therefore: 



N = 





Cpol 



Let us note that this prescription of the noise in the temper- 
ature ILC map added to mitigate the uncertainties due to 
foreground cleaning violates the assumption that the noise 
in temperature is vanishing, used to obtain Eqs. (|2I3I3|) from 
Eq. ©. 

Two masks are considered: one for T and one for Q 
and U. Monopole and dipole have been subtracted from the 
observed ILC map through the HealPix routine remove- 
dipole (jGorski et al.l 120051 '). The same data set has been 
used for the WMAP 7 yr po wer spectrum re - analy sis by 
the QML esti mator BolPol in Gruppuso et alj (|2009 | ) (a nd 
further used in lPaci et all l|201Ch ; iGruppuso et al.l ^Olll D. 



COSMOLOGICAL PARAMETERS 
EXTRACTION 



lLarson et al. (|2010h for the low I analysis. These frequency 
channels have been co-added by inverse noise covariance 



We use CosmoMC (|Lewis and B ridle 2002) in order to com- 
pute the Bayesian probability distribution of model param- 
eters. The pivot scale of the primordial scalar and tensor 
power spec tra was set to fc« = 0.017 Mpc -1 , as recom- 
mended by ICortes. Liddle and Mukherieel (2009). We vary 
the physical baryon density Obh 2 , the physical cold dark 
matter density fl c h 2 , the ratio of the sound horizon to the 
angular diameter distance at decoupling 8, the reionisation 
optical depth t, the amplitude and spectral index of curva- 
ture perturbations ns and log 10 [10 10 A s ]. We assume a flat 
universe, and so the cosmological constant for each model 
is given by the combination J1a = 1 — — ^c- We set the 
CMB temperature T C mb = 2.725 K (|Mather et all Il999h 
and the primordial helium fraction to j/h c = 0.24. We as- 
sume three neutrinos with a negligible mass. In order to 
fit WMAP data, we use the lensed CMB and we follow the 
method implemented in CosmoMC consisting in varying a nui- 
sance parameter Asz which accounts for the unknown am- 
plitude of the thermal SZ contribution to the small-scale 
CMB data points assuming the mo del of . We use a CAMB 
l|Lewis. Challinor and Lasenbvl [20001 ) accuracy setting of 1. 
We sa mple the posteri or using the Metropolis-Hastings algo- 
rithm (|Hastingslll970h at a temperature T = 1, generating 
four parallel chains and imposing a conservative Gelman- 
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Figure 1. Marginalized 68% and 95%-credible contours for 
(r ,n s ) (left panel) and (n s , Hj^f) (right panel) as estimated by 
the WMAP 7 year full likelihood (red lines) and by the BoPix 
plus WMAP 7 year high i likelihood (black lines). 

Rubin convergence criterion (|Gelman and Rubinl I1992D of 
R-l< 0.005. 

With the settings specified above we extract cosmo- 
logical parameters with the latest WMAP likelihood code 
(version v4pl) available at http://lambda.gsfc.nasa.gov/ as 
benchmarks. We prefer to not quote the estimates for the 
cosmological parameters performed by the WMAP team 
since the conventions and the CAMB ve r sion might dif- 
fer fro m those used in lLarson et alj (|2010T l ; iKomatsu et alj 
(|2010h . 

We then extract cosmological parameters by substitut- 
ing the WMAP low-<? likelihood approach with BoPix. In 
doing this we implicitly use the WMAP inputs in polar- 
ization at iV s ide = 16 as described in Section III and not 
those contained in the WMAP likelihood routine publicly 
available. Since temperature and polarization are treated 
at the same resolution by BoPix, we include the WMAP 
high £ likelihood starting at £ = 31 both in temperature 
and temperature-polarization cross-correlation when using 
BoPix, unless otherwise stated. Unless otherwise stated, in 
BoPix we vary the Ce up to £ = 30 and we use the publicly 
available file test_cls_v4.dat as a fiducial power spectrum to 
complete the full covariance at low resolution from £ — 31 
to £ = 64, as done for temperature only by the WMAP pixel 
likelihood. 

We find small differences in the estimate of the cos- 
mological parameters by substituting BoPix to the WMAP 
low-*? likelihood, as reported in Table I. The main differ- 
ence between the estimate of the cosmological parameters 
derived by our alternative low-t? likelihood code and the one 
obtained with the WMAP approach is in the spectral index 
n s : we obtain a value for n s which is 0.86cr lower than the 
WMAP one. This change would lead to quantitative differ- 
ences in the evidence against the Harrison- Zeldovich of the 
WMAP 7 yr data. However, also the other directly sampled 

2 Note that the small differences of our results with the full 
WMAP 7 year likel iho od with respect t o the results reported by 
lLarson et al ] l|2010h or IKomatsu et al . (2010) might be ascribed 
to the different version of RECFAST used, different tools for ex- 
tracting cosmological parameters or different conventions, such as 
the pivot scale fc*. 



cosmological parameters differ from the WMAP estimate in 
about 0.5cr, pointing towards values higher for the physi- 
cal CDM abundance Q, c h 2 and the amplitude of scalar per- 
turbations As and smaller for the baryon physical content 
Qbh 2 and optical depth r. As a derived parameters, we have 
a higher value for the matter content Qm and as, smaller 
for the present Hubble rate Ho. We show more details about 
these different estimates in the two-dimensional plots of Fig. 
[5] These differences seems robust to the change in the multi- 
pole transition to the high likelihood approximation and to 
the change of the fiducial model to complete the covariance 
at low resolution. Special mention should be made for the 
case in which we do not consider It = £p , but we adopt the 
same £ T = 30 and £ P = 23 adopted by the WMAP team, 
but with BoPix for low resolution: the differences with re- 
spect to the estimates by the full WMAP yr likelihood are 
slightly smaller than in the case of It = £p = 30, as can 
be seen in Table I. This means that differences we find are 
not fully due to the different threshold multipoles for po- 
larization adopted in the two low-<? likelihood approaches. 
No appreciable differences are noticed by constructing the 
signal covariance matrix up to 3JV g ide instead up to 4A^id e - 
This can be understood since this different prescription in 
constructing the signal covariance matrix is damped by the 
Gaussian smoothing in intensity and is much below the noise 
in polarization. 

We have performed a further test excluding Asz, just 
for code comparison. We find a smaller discrepancy between 
the estimates for the cosmological parameters and the best- 
fits from the two likelihood approaches when the nuisance 
parameter Asz is omitted (i.e. fixed to zero). This addi- 
tional foreground parameter Asz is not well constrained by 
WMAP, but it contributes to the shape of the final likelihood 
and to the marginalized values of the parameters (shifting 
slightly the value of n s , for instance). 

Most of these small differences reported in the estimate 
of the cosmological parameter s interfere destructively be - 
cause of the cosmic confusion (|Efstathiou and Bondll998l ) 
and the best-fits Ce from the two likelihood analysis agree 
very well. We present the CMB bestfit Cj T obtained by 
BoPix in combination with the WMAP 7 high-<? likelihood 
in comparison with those obtained by the full WMAP 7 like- 
lihood in Fig. [S] The difference in the best-fit Ce is consistent 
with the different central values for the cosmological param- 
eters displayed in Table [5] We have checked that the best-fit 
C t obtained in this work by the full WMAP 7 likelihood 
has A(— 2 log £wmap) = —7.42 with respect to the refer- 
ence WMAP 7 test_cJs_v4.dat; the best-fit Ce obtained in 
this work by BoPix in combination with the high-*? WMAP 
7 likelihood provides a better fit, with A(— 2 log £wmap) = 
—7.75 with respect to the reference WMAP 7 test-ds-v4.dat. 

Since the pixel approach and the Gibbs approach con- 
tained in the WMAP likelihood package give quite consistent 
results, the difference we find by using BoPix is more likely 
due to the asymmetric treatment of T and (Q, U) adopted 
in WMAP low-<? likelihood. We have therefore tested BoPix 
against the WMAP likelihood within the same range of mul- 
tipole, i.e. up to £ — 30: BoPix has been run on the low- 
resolution WMAP 7 yr N silic = 16 products varying Cj T , 
Cf E , Cj E up to £ = 30 and compared to the likelihood ob- 
tained by the WMAP 7 yr pixel based routine plus the high-*? 
likelihood value for TE from £ = 24 to £ = 30. In this way 
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Parameter 



WMAP 7 
likelihood 



WMAP 7 
likelihood 



e T = t P = 30 i T = i P = 24 e T 



36 



Ip 



30 

: 23 



Different 
Fiducial 



(Pixel) 



(Gibbs) 



100 Cl b h 2 

n c h 2 



log[10 10 A 



2.250 ±0.056 

1H4+0 0054 
u - lxl -0.0053 

0.089 ±0.015 

n 11fi +0.033 
°- llu -0.032 



2.252 H 



f0.057 
'-0.056 

0.1114 ±0.0055 



o.o89i2;2^ 

o.969l2;2i 3 4 

3.116 ±0.033 



2.213 ±0.055 

1145+ - 0055 
u.nio_ 0056 

085+ ' 015 
u.uso_ 014 

0.956 ± 0.014 

3.130 ± 0.033 



2.215 ±0.055 
0.1142 ±0.0055 

o.ossig^ 
0.95712:2^ 

3.128 ±0.032 



2.224 



+ 0.057 



0.058 
+0.0055 
0.0056 
+0.014 
0.015 
Q54+0- 013 

U.9D1_ Q14 
n , OO+0.032 

o.ioo_ 033 



0.1152 
0.085 



2.224 ±0.055 

n i 19 o+0.0058 
u.iizo_ 0057 

0.087 ±0.015 

0.959 ±0.014 

3.122 ±0.033 



2.212 ±0.058 
0.1144 ±0.0056 
0.085 ±0.015 
956+ ' 014 

U.9OD_ 013 

3.128 ±0.033 



Ho 

(78 



270+ - 027 
U - Z ' U -0.028 

70.71^ 
0.811 ±0.029 



0.269 ±0.028 

7O.7I24 
0.811 ±0.029 



0.289 



+0.031 
-0.030 

68.9I2.4 

0.820 ±0.029 



n 9HS+0-030 
0.288_ Q31 

69.1 ± 2.4 

n siq+O- 028 
u.oi3_ 029 



0.294 ±0.031 
68.5 ± 2.4 



0.822 



+0.030 



n 977 +0.031 
u ' z ' ' -0.030 

69.9 ±2.5 
n »in+ 031 

U.810_ ao30 



0.289 ±0.030 

68.913-4 

n R9n+ 029 
u.8Ai_ a030 



Table 1. Mean parameter values and bounds of the central 68%-credible intervals for the cosmological parameters estimated by the 
WMAP 7 year full likelihood (second and third column) and by the BoPix plus WMAP 7 year high I likelihood for different transition 
multipolcs It = ^P (fourth, fifth and sixth column), for It 7^ t-P and different fiducial theoretical power spectrum to complete the signal 
covariance matrix in BoPix (last column). Below the thick line analogous mean values and bounds are presented for derived parameters. 



, f 3000 r 






Figure 2. Comparison of the best-fit Ci obtained by BoPix in 
combination with the WMAP 7 high-£ likelihood (solid) vs. the 
WMAP 7 full likelihood (dashed) is shown in the top panel. To 
make the difference more visible, the relative difference between 
the C( bestfits in temperature is shown in the bottom panel. Note 
that the differences are well within the cosmic variance. 



we subtract the same high-£ likelihood information from hy- 
brid runs presented in Table I. By assuming n b h 2 = 0.02246, 
Q, c h 2 = 0.1117 and sound horizon 6 = 1.03965, we obtain 
results quite consistent with the hybrid ones: a slight smaller 
value in the estimate of r and ns and a larger one for As- 



Figure 3. Marginalized one-dimensional probabilities for r, n s 
and logllO 10 ^] as estimated by the WMAP 7 year full likelihood 
(red lines) and by the BoPix plus WMAP 7 year high I likelihood 
(black lines). See text for further details. 



6 OTHER EXTENDED COSMOLOGICAL 
MODELS 

We now consider few cosmological models beyond the 
ACDM model which can be constrained by WMAP 7 year 
data only. We consider only the baseline /trans = 30 and 
all the other settings consistently with the previous section, 
unless otherwise stated. 

Gravitational Waves. 
We consider all inflationary models which can be described 
by the primordial perturbation parameters consisting of the 
scalar amplitude and spectral index {As, ns), and the tensor- 
to-scalar ratio r. In canonical single-field inflation, in the 
slow-roll limit, the tensor spectrum shape is not indepen- 
dent of the scalar one. We will consider a tensor spectrum 
with a tilt ut = — r/8, as predicted for canonical single-field 
inflation at first-order in slow-roll. 

Our marginalised 68%-credible interval for the scalar 
spectral index is given by ns = 0.977l2'o2i (to be compared 
with the result of the previous section, ns = 0.955 ± 0.014), 
half a sigma redder than the result we obtain by the full 
WMAP 7 year likelihood 0.987 ± 0.020. 
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0.1073 ± 0.0063 

091+ - 015 
u.u»i_ Q 014 

0.987 ±0.020 
3.093 ± 0.038 


2.270 ± 0.073 

0.10991^066 
0.0871°:* 

g77 +0.020 
u - yl ' -0.021 

3.102 ±0.039 


r 


< 0.34 


< 0.36 


Qm 
H 


246+ - 031 
u - z ™-0.032 

73.2 ± 3.2 


262+ 035 
u.^o^_ 036 


0"8 


0.797 ± 0.033 


0.805 ± 0.033 



Table 2. Mean parameter values and bounds of the central 68%- 
credible intervals for the cosmological parameters including the 
tensor-to-scalar ratio estimated by the WMAP 7 year full likeli- 
hood (left column) and by the BoPix plus WMAP 7 year high £ 
likelihood (right column). For the tensor-to-scalar ratio r the 95%- 
credible upper bound is quoted. Below the thick line analogous 
mean values and bounds are presented for derived parameters. 
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Figure 4. Marginalized 68% 95% contours for (n a , r) as esti- 
mated by the WMAP 7 year full likelihood (dashed lines) and by 
the BoPix plus WMAP 7 year high I likelihood (solid lines). The- 
oretical predictions of few popular inflationary models (including 
reheating uncertainties where appropriate) are displayed. 

At 95% confidence level, our result for the tensor-to- 
scalar ratio is r < 0.36, fully consistent with the result we 
obtain from the full WMAP 7 year likelihood, i.e. r < 0.34. 
Let us note that, differently from the WMAP low-^ likeli- 
hood code, the value obtained by BoPix include BB polariza- 
tion in the construction of the covariance at low resolution. 
The differences in the (ns,r) are shown in Fig. 2 and are 
mainly due to a shift of the constraints at smaller values 
for ns, as occurs for the standard ACDM model discussed 
in the previous section. Theoretical predictions of few pop- 
ular inflationary models (including reheating uncertainties 
where appropriate) are displayed. One of the phenomeno- 
logical differences from the different constraints would be a 
minor tension for a massless self-int eracting inflaton mode l 
with WMAP 7 ye ar data only (see iKomatsu et all (|2010T l; 
iFinelli et al] (|2010T l as examples for an higher tension of the 
X4> 4 potential with observations when additional cosmolog- 
ical data sets are added to WMAP). 

Runnvng of the scalar spectral index. 
In this subsection we consider the variation of the scalar 
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Table 3. Mean parameter values and bounds of the central 
68%-credible intervals for the cosmological parameters including 
the running of the scalar spectral index n run estimated by the 
WMAP 7 year full likelihood (left column) and by the BoPix plus 
WMAP 7 year high i likelihood (right column). For the running 
of the scalar spectral index n run the 95%-credible upper bound 
is quoted. 
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Figure 5. Marginalized 68% and 95%-credible contours for 
(n s ,n r un) as estimated by the WMAP 7 year full likelihood (red 
lines) and by the BoPix plus WMAP 7 year high I likelihood 
(black lines). 

spectral index with wavelength, i.e. we allow n run to vary in 
the range [—0.2, 0.2]. Our marginalised 95%-credible interval 
for the scalar spectral index is given by —0.065 < n run < 
0.042, which has to be compared with the result we obtain 
by the full WMAP 7 year likelihood -0.074 < n run < 0.030. 
The results are fully consistent between each other and with 
the hypothesis of no wavelength dependence of the scalar 
spectral index. 

Neutrino Mass. 

In this subsection we constrain the total mass of neu- 
trinos m v = 94Q„h 2 eV, allowing to vary the fraction of 
massive neutrino energy density relative to the total dark 
matter one /„ = Q^/Qum- At 95% confidence level, our 
result for the fraction of massive neutrinos is /„ < 0.113, 
whereas we obtain f„ < 0.094 from the full WMAP 7 year 
likelihood. The resulting neutrino mass bound at 95% confi- 
dence level is ^ m v < 1.4 eV, compared to 1.1 eV obtained 
from the full WMAP 7 year likelihood. 

Cosmological Birefringence. 
Since one of the main differences between the WMAP low 



© 2010 RAS, MNRAS 000, [TV?? 



Cosmological Parameters from a re-analysis of the WMAP-7 low resolution maps 7 



PaxEiiiictcr 


WMAP 7 


BoPix plus 




likelihood 


WMAP 7 high I likelihood 


100 n b h 2 
n c h 2 


9 91Q +0.062 
z.zia_ 060 

~ -, „„ 4-0 0071 
o.ll 1 1 _ .0073 


2.174 ±0.061 

„ ^ ~~^4-0 0081 
u,i " u -0 .0080 


T 


n no^+0.014 

°- 087 l0.015 


0.082 ± 0.014 


n s 


0.960 ± 0.016 




log[10 10 A s ] 


3.120 ±0.032 


3.134 ±0.033 


fu 


< 0.094 


< 0.113 


Qm 
H 
era 


, 2g +0.057 
65.7ilj 

n 71 9 +0.073 
U - ' 1 -0.74 


u -°' -0.072 

62.8+t? 
69 J +0 b87 

U.D9O_ 0g3 




< 1.1 cV 


< 1.4 eV 



Table 4. Mean parameter values and bounds of the central 68%- 
credible intervals for the cosmological parameters including the 
total mass of the neutrinos estimated by the WMAP 7 year full 
likelihood (left column) and by the BoPix plus WMAP 7 year high 
I likelihood (right column). For the total mass of the neutrinos 
m„ the 95%-credible upper bound is quoted. 
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Figure 6. Marginalized 68% and 95%-credible contours for 
Ei/" 1 * i^M^ 2 ) (left panel) and (n a ^22, v m v) (right panel) as 
estimated by the WMAP 7 year full likelihood (red lines) and by 
the BoPix plus WMAP 7 year high I likelihood (black lines). 

resolution likelihood code and BoPix is the treatment of the 
polarization sector, we now wish to analyze an extended 
cosmological model different from ACDM only in (Q, U) 
and the relative cross-correlation with the temperature. Cos- 
mological birefringence refers to a non-vanishing interac- 
tion oc cf>F l u,F fl1 ' between photon and a cosmological evolv- 
ing pseudo-scalar <f), which would generate non-vanishing 
TB and EB correlations (|Lue. Wang and Kamionkowsk.il 
Il999l ) through a rotation a of the polarization plane of 
CMB photons along their path from the last scattering sur- 
face to the observer. The resulting polarization and cross 
temperature-polarization spectra would encode the partic- 
ular redshift dependence of the parity violation interaction 
iFinelli and Galavernil l|2009h . However, a phenomenological 
shortcut exists, commonly used in the literature and also 
adopted by the WMAP team, and consists to neglect the 
redshi ft dependence of a and simple predict the power spec- 
tra as lLue. Wang and Kamionkowskil (|l99Qt ): 

^EE,obs r^EE 2 /n \ 

C t = C t COS (2a) , 

C e = C e sin (2a) , 



,~tEB,obs 1 y~tEE • (a \ in\ 

C e = -C e sm(4a) , (7) 

nTE,obs nTE \ 

C e = Ci cos(2a) , 

^TB,obs syTE ■ /r, \ 

C e = C e sm(2a) . 

The above formulae are valid when the primordial B-mode 
polarization is negligible, which is assumed in this paper. 

We have therefore sampled a in radiants with a flat 
prior [—0.5 0.5] plus the other six cosmological parameters 
of the ACDM model by inserting Eqs. (J7J). Our marginalised 
68% (95%)-credible interval for a is a = -1 °-3+°,°; 7 1 2 ,"'! 
in agreement with the full WMAP 7 year likelihood result 
which we find a — — 1 ".O^Qo'gl^o 3 Either the result us- 
ing BoPix or the one based on the full WMAP 7 year like- 
lihood are consistent with vanishing cosmological birefrin- 
gence at 95% CL just by assuming the statistical uncer- 
tainty, and the agreement increases by using the systematic 
uncertainty, wh i ch is e stimated as 1 ° .4 by the WMAP team 
iKomatsu etaTI (|2010l ). 

Since the weight of the high-^ TB likelihood plays a 
relevant role in these constraints we have also considered 
the case in which this is not taken into account. Such 
setting which emphasizes the role of polarization on large 
angular scales would be relevant to show clearly the po- 
tential differences between BoPix and the WMAP pixel 
likelihood code. On using only low resolution products to 
constrain cosmological birefringence, by using BoPix on 
A^idc = 16 resolution Q , U maps and matrices we obtain 
a = -4°.2^j„'j^ 7 o 5 2 , still in agreement with the values we 
find by the WMAP 7 likelihood on N sidc = 8 resolution Q , U 

_[_ Q O g I 1 A O Q 

maps and matrices a = —0 .2_ 30 ' 6 _ go g Although with 
larger uncertainties, our results agree with vanishing cosmo- 
logical birefringence at 95% CL, without invoking system- 
atic uncertainties. Note also that o ur result agrees with th e 
analysis on large angular scales by iGruppuso et alj l|2012h , 
which has much tighter constraints probably because all the 
cosmological parameters except a are kept fixed. 

The full posterior likelihood for a and its two dimen- 
sional contour in combination with the optical depth r are 
shown in Fig. (JSJ). Fig. (JS| also shows that no degeneracy 
between r and a is observed in WMAP 7 yr data. Note that 
the slight preference at 68% CL for negative values of a 
when using only BoPix on low resolution products is consis- 
tent with the WMAP 7 yr TB and EB power spectra QML 
estim ates at i < 30 and presented in lGruppuso et all (1201 ll . 

En!). 



7 CONCLUSIONS 

We have performed an alternative estimate of the cosmo- 
logical parameters from WMAP 7 year public data, by sub- 
stituting the WMAP 7 low-^ likelihood with a pixel like- 
lihood code which treats (T,Q,U) at the same resolution 
without any approximation. We have used this code at the 
HEALPIX resolution N^c = 16 on foreground cleaned pub- 
lic data, therefore increasing the resolution of the pixel based 
polarization products used in our extraction of the cosmo- 
logical parameters with respect to the WMAP standard one. 
We have consistently increased the transition multipole from 
t = 24 to £ = 31 for the high-^ WMAP 7 year temperature- 
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Table 5. Mean parameter values and bounds of the central 68%- 
credible intervals for the cosmological parameters allowing for an 
effective treatment of cosmological birefringence estimated by the 
WMAP 7 year full likelihood (left column) and by the BoPix plus 
WMAP 7 year high £ likelihood (right column). For the angle a 
defined in Eq. (Q the 95%-credible upper bound is quoted. 
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Figure 7. Marginalized posterior probability for a (left panel) 
and marginalized 68% and 95%-credible contours for (t , a) (right 
panel) as estimated by the WMAP 7 year full likelihood (dashed 
red lines) and by the BoPix plus WMAP 7 year high I likelihood 
(solid black lines). The additional dot-dashed blue line and short- 
dashed pink lines are for the constraints on a from large angular 
scales only obtained by the WMAP 7 year pixel likelihood code 
and BoPix, respectively. 

polarization cross-correlation likelihood and included the 
marginalization over the nuisance parameter Asz- 

With this setting we have found estimates for the cos- 
mological parameters consistent with those obtained by the 
full WMAP 7 year likelihood package, although for some 
parameters the differences are of half a or more. These dif- 
ferences between the two \ow-£ likelihood treatments we 
find are larger than the WMAP 7 yr likel ihood uncer- 
taintie s from tests on simulations reported in lLarson et al.l 
(2010); however, we need to keep in mind that our differ- 
ences between two likelihood treatments are reported for real 
data, with WMAP 7 year beam/points source corrections 
and marginalization on Asz taken fully into a ccount, differ- 
ently f rom the simulation analysis performed in lLarson et al.l 
(2010). It is also important to stress that the WMAP as- 
sumption to have vanishing noise in temperature is explicitly 
violated by the ILC WMAP 7 yr temperature map digested 
by the low-^ likelihood code (the likelihood test on simu- 



lations reported in lLarson et all (|2010h might also miss the 
additon of the 1/jK noise per pixel to the low-resolution tem- 
perature maps). Another important aspect is that the low 
resolution polarization data digested by the two likelihood 
are different since are at different resolution. The difference 
between the two best-fit Ce for ACDM found by the two 
alternative likelihood treatments show a maximum of 4% 
around at £ ~ 10 and oscillate with an amplitude below 1% 
for i > 1000. 

On restricting to the ACDM model the most important 
difference is for the scalar spectral index ns, which decrease 
to 0.956 from the value 0.968 we obtain with the full WMAP 
7 yr likelihood code, i.e. a decrease of 0.86 a. This differ- 
ent value for ns would increase the evidence against the 
Harrison- Zeldovich spectrum from WMAP 7 yr data. This 
difference for ns is consistent with the one between the two 
best-fit Ce and depend only partially from the threshold 
multipole from which the high-i' TE likelihood starts. Other 
previous alternative likelihood treatments also reported the 
most important dis c repancy for the sca lar spectral index 
( Eriks en et al ] |2007l : iRudiord et aljlioosl ). A smaller value 
for ns with respect to the estimate by the full WMAP 7 
year likelihood code, always within 1 <r, is then seen in all 
the extension of ACDM considered here. No major changes 
are found for the 95 % credible intervals for the tensor to 
scalar ratio and for the running of the scalar spectral index. 
A slight degradation has been found for the 95 % credi- 
ble interval on the neutrino mass. The case of cosmological 
birefringence has been taken as a sensitive test for the two 
alternative likelihoods, whose most relevant difference is the 
treatment of polarization on large scales. As expected, a 
slight difference on the posterior of the polarization angle 
a has been found when only low resolution data are used, 
whereas the results are fully consistent when the high-^ TB 
data are added to both likelihoods. 
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